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Abstract 

We study the interference effects between three weakly hnked trapped Bose- 
Einstein condensates (BEC) as a generahzation of the two-component conden- 
sates. Three coupled Gross-Pitaevskii equations (GPE) are used to describe 
the dynamics of the system. The nonsinusoidal oscillation is found as a gener- 
alization of the Josephson effect in superconductivity. The self-trapped effects 
are also predicted in three coupled BEC. Moreover, in general case, the phase 
diagrams of the system are closed only for some special parameters, which 
can be used to determine the interaction parameters between atoms in BEC. 
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Recently, Bose-Einstein Condensate (BEC) has been realized in a dilute and ultracold gas 
of trapped alkali-metal atomsliHl. It opens a new field for studying the physics of BEC, which 
is hoped to improve our understanding on the fundamental concept of quantum mechanics 
and quantum many body problem&i. From the ensuing theoretical and experimental works, 
it is clear that the condensed Bose gas displays the crucial properties of a confined coherent 
matteJi. 

Since the current understanding of BEC is mainly based on the mean-field approximation, 
where a macroscopic wave function is introduced as the order parameter, the study of this 
feature by investigating the interference phenomena should be of great importance&^l. Such 
a description using a macroscopic wave function with a definite phase implies U(l) gauge 
symmetry broken. 

There are many methods suggested theoretically for detecting the interference effects 
between two components of condensates. The one is the investigation using the continuous 
measurement theoryi. It shows that an interference pattern between two condensates can 
be built up dynamically in a single run of an experiment, even though no phases have even 
been assumedB^. Another approach is that, one can set up a potential barrier between 
two trapped condensates by using a far off-resonant laser. If one lowers the potential barrier 
produced by the laser, the atoms would tunnel through the barrier due to weak coupling 
between these two condensates. The interference is developed, somewhat like Josephson 
effect in super conduct ivityEl'0. Similarly, if one switches off the laser potential traps, the 
atoms would expanded ballistically due to the sudden disappearance of the barrier. The 
interference fringes can be observedllUli. It is also interesting experimentally that the inter- 
ference effects can also be tested by means of different hyperfine atomic states in a single 
trap0. 

Among these methods mentioned above, some have been realized in experiments^H^l. In 
Ref. [20], the authors reported that fringes were observed in the density of two overlapping 
condensates, released after switching off the traps. In Ref. [21], the authors measured the 
relative phase of two condensates in different hyperfine atomic states occupying a single trap. 



using Ramsey's method of separated oscillating fieldH. From experimental dataliatJ, the 
phase locking indeed occurs for small separation between condensates, implying the broken 
gauge symmetry. 

The interference effects of two component condensates have been studied in a variety of 
cases, but there are very few on threei. In this paper, we would discuss the interference 
effects between three coupled BEC. 

Any generalization on this topic would increases the difficulty significantly. The reward is 
that the more complex nonlinearity and the more parameters in the basic equations would 
bring richer dynamics and more compound structure into the model. Experimentally, a 
far off-resonant laser barrier can be used to split a trapped condensate into two, three or 
even more well-trapped condensates. Lowering the laser intensity allows atoms tunneling 
through the barrier. Then interference will be built up between these condensates. The 
optical confinement method also has the possibility to create a three component BEC in 
different hyperfine atomic states^. 

Due to the U(l) gauge symmetry broken, it is crucial to introduce a global phase to the 
order parameter, the macroscopic wave function \E'(r, t) to describe Then, the 

wave function '^(r,t) at T = can be factorized aJii@ii 

^(r,t)=^(t)$(r), (1) 

which satisfies a nonlinear Shrodinger equation, the Gross-Pitaevskii equation (GPE)0, 

= _|!.v2vl/(r,t) + [Vtrapir)+go \^ (r , t)f]^ (v , t) , (2) 

where Vjrap(i") is the external trap potential and go = Anh'^a/'m^ m the atomic mass and a 
the scattering length of the atom-atom interaction. When we study the interference of three 
weakly coupled BEC in traps 1,2 and 3, the dynamics of the system is described by three 
GPE, coupled with tunneling matrix elements. It can be written as a three-state modefi 

ih^ = + UiN,)ij, - Kij2, (3) 
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= {El + ^72iV2)^2 - i^^i - i^^3, (4) 

ih^ = {El + U^N,)^^ - K^2, (5) 



with 



"01,2,3 = ^A^i,2,3 exp(i6'i,2,3(t)), (6) 

where iVi 2,3 and 6'i^2,3 are the atom numbers and (time-dependent) phases for BEC in trap 
1, 2 and 3, respectively. K is the couphng matrix element between trap 1 and 2 or trap 
2 and 3 (for simplicity, we only consider the case, that these two coupling matrix elements 
are the same and there is no coupling between 1 and 3). -E123 are the zero-point energies 
in each well, f/i, 2,3^1,2,3 are proportional to the atomic self-interaction energies, induced by 
the overlap of the spatial part of the macroscopic wave function, and the total number of 
atoms Nt = Ni + N2 + is a. constant. 

From Eqs. (|), @, (^) and (||), we get four independent equations in terms of the phase 
difference and the atom population imbalance 

= {El - El) + {U,N, - U2N2) + COS01 + i^y^cos(0i - <f>2), (7) 

= {El - El) + (f/iiVi - U3N3) - ir^cos</>i + ir^cos(0i - <P2), (8) 

hd{Ni-N2) 



2 dt 
hd{Ni-N3) 



2KJNiN2sm<f)i - K J N2N3 sm{<f)i - ^2), (9) 



KJNiN2Sin(P^ + KdN2N3sm{(Pi - ^2), (10) 



2 dt 

where 0i = 6*2 — Oi and 02 = ^3 — ^i- 

In order to get an understanding on the novel phenomena of interference effects between 
three trapped BEC, let us consider an ideal case: the totally symmetric case, in which, 
Ui = U2 = U3 = U, El = El = El, and the atom population imbalances of traps 1, 2 and 
of traps 2, 3 are the same, that is Ni — N2 = — N2, or Ni = N^. Then, we can simplify 
Eqs. (0), (i, (i and (0) as 



dd^ N-i - A^o /iVT 

= [/(TV, - A^,) + K^=^cos4., + Ay^cos^,. (11) 

hd{Ni-N2) 



-3KJNiN2sm(j)i. (12) 
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After rescaling the time variable and the population imbalance, 2Kt/h —* t, z = ^ jv^^^ ; 
Eqs. ( p!TD and ([T2|) can be written as 

dz 



— = -Vl -z-2z^ sin 0, (13) 
^ = AzH =cos(^, (14) 

where A = = 0i. 

If we treat z{t) and 0(t) as a pair of canonically conjugate variables, it is not difficult to 
construct a Hamiltonian 

A 



H = -z"^ -Vl-z-2z^cos(j), (15) 



satisfying the canonical equations, 
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In Eq. (^), 2; < 0.5, consistent with the totally symmetric conditions. The fluctuations of 
the phases and the atom numbers must be small and has not been included herelll. 
Now, let us discuss Eqs. ([T3|) and (|14D in more details. 

(1) zero-phase mode. With initial conditions z{0) = 0.3 and 0(0) = 0, we can solve Eqs. 
(p!3|) and (|l^) numerically in zero-phase mode. Some results are given in Fig. |l]. A takes on 
the values 10, 20, 38, 38.25 and 39 for Fig. 1(a), 1(b) ,l(c), 1(d) and 1(e), respectively. The 
anharmonic nonsinusoidal oscillations are the generalized sinusoidal Josephson oscillationsii. 
z is asymmetric about its zero point because of the existence of the third trap, {z{t)) 7^ 0. 
As A increases, the oscillations around z = become anharmonic [Figs. 1(a), 1(b) and 1(c)]. 
When A exceeds a critical value Ac = 38.25, the populations oscillates around some nonzero 



values, and z{t) > [Fig. 1(e)]. This population imbalance self-trapped effect is somewhat 
like the pendulum bob swing over the (p = tt vertical orientation for sufficiently large initial 
values. This state can be achieved from different approaches. In zero-phase mode, there is 
a pair of eigenf unctions of GPE, z = and (p = it, and the ground state energy is E = 1. If 
the self-trapped effect takes place, z when (p = n, then 

A 



Hiz{0),m) = ^jm' - VI - ^(0) - 2^(0)2 cos 0(0) > 1, (18) 
so the critical values depend on the following condition 



1 + Jl - z(0) - 2^(0)2 cos 0(0) 

Here, ^(0) = 0.3 and 0(0) = 0, so = 38.25. 

(2) 7T-phase mode. If we select the initial conditions as 0(0) = vr and z{0) = 0.3, we 
can discuss the interference dynamics of vr-phase oscillations from Eqs. (p!3| ) and (p!^. The 
numerical results are shown in Fig. 0. It is clear that the nonsinusoidal oscillations in vr-phase 
mode are similar to that in zero-phase mode [Fig. 2(a)], but the self-trapped effects are not. 
There are two kinds of self-trapped effects in vr-phase mode {z(t)) < z{0) and {z{t)) > z{0) 
as given in Figs. 2(d) and 2(f). According to the degenerate GPE eigenstates that break the 
z symmetry, we can get the critical value of Ac for these two different self-trapped effects 

A. ^ " ' . (20) 

2^(0)^1-^(0) -2^(0)2 

here, z{0) = 0.3, so A^ = 5.08. If A exceeds the value of A^,, the system goes from the first 
kind of self-trapped state [Fig. 2(d)] into the second kind of self-trapped state [Fig. 2(f)]. 

(3) z-(j) phase diagrams, z and are two canonically conjugate variables in Eqs. (p!3D 
and (p!^, their dynamical behavior in z-(j) phase diagrams are shown in Fig. ^. In Figure 
3(a), A = 10 and 0(0) = 0, it is the zero-phase mode. For z{0) = 0.20, and 0.40, the 
trajectories are closed, asymmetric curves, corresponding to the generalized nonsinusoidal 
oscillations. For z{0) = 0.487, it is the critical case (here, for the zero-phase mode, we 
should solve the critical value of z from Eq. ([T9|) with A = 10). And for z{0) = 0.498, the 



running self-trapped case, z{t) is locked but (j)(t) is unlocked, — oo < (j)(t) < oo. In Fig. 3(b), 
A = 7 and 0(0) = tt, the vr-phase mode, the system is self-trapped for all values of z{0). For 
z{0) = 0.08, it is the running-mode self-trapped case, is unbounded, — oo < 0(t) < oo, 
while for larger z{0), above a critical value of Zc = 0.111 according to Eq. ( [20| ) (here, for 
the TT-phase mode, we should solve the critical value of z from Eq. ( PP] ) with A = 7), 0(t) 
is also localized around it. So in Fig. 3(b), there are two kinds of self-trapped effects for 
z{0) = 0.08 and for ^(0) = 0.111, 0.30 and 0.37 in vr-phase mode, respectively. 

The self-trapped effects in coupled BEC are the results of the nonlinearity of interatomic 
interaction in GPE and the long-range quantum coherence of a macroscopic number of 
atoms. 

(4) z-(f) phase diagrams, the general case. From Eqs. (^, (||), (|) and (p!OD generally, 
after rescaling the atom population imbalance variables zi = '^'^j^^^ and Z2 = ^ V/^^ ' 
time 2Kt/h — > t, we can get four coupled dynamical equations 



1 3zi 1 /zi + 1 - 22:2 

: cos 01 + - W ; COSf 01 - 02 J 



dt 2 ^(^i + 22 + l)(^2 + l-2zi) 2\ Z2 + 1-2Z, 



+ ^A{zi + Z2 + l)-^A{z2 + l-2zi), (21) 

902 1 IZ2 + 1- 2Zi 1 IZ2 + 1-2ZI , , , , 

1^ = ~oV ^ — COS01 + -\ -r— —COS 01 -02 

Ot 2 \l Zi + Z2 + 1 2 V 2^1 + 1 — 2^2 

+ ^A{z^ + Z2 + l)-^A{z, + l-2z2), (22) 
dzi 2 



^^izi + Z2 + l){z2 + l-2zi)sm(f)i 
-1\/{Z2 + 1 - 2z,){z^ + 1 - 2^2) sin(0i - 02), (23) 
^ = -^'\/{zi + Z2 + l)(z2 + 1 - 2zi)sin0i 



+ ^^{Z2 + 1- 2zi){zi + 1 - 2^2) sin(0i - 02). (24) 

Here, for simplicity, symmetric parameters are also considered, Ui = U2 = U3 = U, = 
E2 = -E3, and A = Besides the nonsinusoidal oscillations and self-trapped effects, 

we can find another novel phenomenon between three trapped BEC. With some initial 
conditions (0i(O), 02(0), -2i(0), ,22(0) and A), the trajectory of z-cp phase-diagrams (-2i-0i, 



Z2-<p2) are not closed. They are unclosed for general initial conditions, but closed for some 
special initial values. As shown in Fig. ^, we fix the values of 0i(O) = 0, 02(0) = 0, 
^;^(0) = 0.3 and ^2(0) = 0.1, and let A varied. When A = 2 [Fig. 4(a) and 4(b)] and 
A = 9 [Fig. 4(e) and 4(f)], the z-cj) phase-diagrams are unclosed. But for A = 6.8 [Fig. 
4(c) and 4(d)], they are convergent and closed. These phase-diagrams exhibit a transition 
from unclosure to closure, then to unclosure again during changing the values of A. Only a 
special value of A can make phase- diagrams close under some values of 0i(O), 03(0), -2i(0) 
and -^2(0). This property can be used to determine the value of A in experiments. 

This feature of phase-diagrams originates from two aspects. One is the nonlinearity of 
the interatomic interaction. It is known that the phase-diagram of a nonlinearity system is 
always not closed0. Eqs. (pi]), (^2]), ( ^3]) and (^4]) are four coupled nonlinear equations, 
so the phase-diagrams are generally unclosed. The other is the interaction between three 
trapped BEG. The interference of BEC in a macroscopic number of atoms determine the 
properties of phase-diagrams. The detail study on this complicated problem would be a 
special topic elsewherelil. 

In summary, the interference of three trapped BEC induces population nonsinusoidal os- 
cillations which is a generalization of sinusoidal Josephson effect in superconductivity. The 
quantum self-trapped effects occur in three BEC for both zero-phase mode (one kind of 
self-trapped state) and vr-phase mode (two kinds of self-trapped states) when parameters 
exceed some critical values. Any observation of the predicted interference effects would un- 
ambiguously prove the existence of gauge symmetry broken and the relative phase between 
condensates. This consideration, concerning the interference between three Bose conden- 
sates, allows the numerical simulation of some recent experiment sil. From calculations, we 
have inferred that the convergence of the phase diagrams is sensitive to the adjustment of the 
system parameters. The phase-diagrams are not closed in general cases, but for some special 
conditions, the phase-diagrams are convergent and closed, which can be used to determine 
the interaction parameter between atoms of BEC in experiments. 
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FIGURES 

FIG. 1. Rescaled population imbalance z{t) versus dimensionless time variable t, with initial 

conditions z(0) = 0.3 and 0(0) = in zero-pliase mode. The interaction parameter A takes the 
values (a) 10, (b) 20, (c) 38, (d) 38.25 and (e) 39. 

FIG. 2. Rescaled population imbalance z{t) versus dimensionless time variable t, with initial 
conditions 2;(0) = 0.3 and (f)(0) = tt in 7r-phase mode. The interaction parameter A takes the values 
(a) 2, (b) 4, (c) 5, (d) 5.05 (e) 5.08 and (f) 5.1. 

FIG. 3. Phase diagrams of population imbalance z versus phase difference 0. (a) and (b) are 
corresponding to zero-phase mode <p{0) = with A = 10 and 7r-phase mode <p{0) = it with A = 7, 
respectively. The values of z{0) are as given. 

FIG. 4. Phase diagrams of population imbalance versus phase difference in general case, with 
0i(O) = 0, (/>2(0) = 0, zi{0) = 0.3 and ^2(0) = 0.1. (a), (c) and (e) for zi — (pi phase diagrams, 
while (b), (d) and (f) for Z2 — ^2 phase diagrams. The dimensionless time —65 < i < 65. A takes 
the values 2 for (a) and (b), 6.8 for (c) and (d), 9 for (e) and (f). 
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